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CONVERSION  FACTORS,  NON- SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  can  be  converted  to  SI  (metric)  units  as  follows: 


_ Multiply _  By 

inches  25.  A 

feet  0 . 3048 

pounds  (force)  4.448222 

pounds  (mass)  per  cubic  foot  16.01846 

Btu  1.055 

Fahrenheit  degrees  5/9 


_ To  Obtain _ 

millimetres 

metres 

newtons 

kilograms  per  cubic  metre 
kiloj  oules 

Celsius  degrees'*  or  kelvins 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  -  (5/9) (F  -  32).  To  obtain  Kelvin  readings, 
use  K  -  (5/9) (F  -  32)  +273.15. 
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PREDICTING  CONCRETE  SERVICE  LIFE  IN  CASES  OF 


DETERIORATION  DUE  TO  FREEZING  AND  THAWING 


PART  I :  INTRODUCTION 

1.  The  Coras  of  Engineers  (CE) ,  in  support  of  its  navigation  and  flood 
control  missions,  must  maintain  the  integrity  of  a  large  number  of  mass  con¬ 
crete  structures.  It  often  is  necessary  to  evaluate  the  condition  of  the 
concrete  in  these  structures  to  determine  the  magnitude  and  significance  of 
specific  problems,  to  assess  the  need  for  remedial  action,  or  even  to  provide 
a  basis  for  no  action  in  conditions  that  might  seem  to  indicate  some  action. 

A  consistent  and  uniform  procedure  for  predicting  the  service  life  of  these 
concrete  structures  is  necessary  to  properly  allocate  public  funds  for  solu¬ 
tion  of  identified  problems .  Such  a  procedure  should  reflect  the  current 
understanding  of  mechanisms  which  cause  premature  deterioration  in  concrete. 

It  should  incorporate  available  information  concerning  the  concrete,  environ¬ 
mental  factors,  and  the  current  state  of  deterioration.  In  addition  to  incor¬ 
porating  what  is  known  about  these  factors,  the  prediction  of  service  life 
should  reflect  the  uncertainties  in  knowledge  of  the  foregoing  factors. 

2.  A  procedure  for  predicting  the  service  life  of  concrete  subject  to 
damage  by  freezing  and  thawing  which  meets  the  above  criteria  has  been  devel¬ 
oped  and  reported  previously  (Bryant  and  Mlakar  1988) .  This  procedure 
addresses  both  the  known  and  the  uncertain  qualities  of  the  relevant  material 
properties,  environmental  factors,  and  model  of  degradation  due  to  freezing 
and  thawing  using  a  probabilistic  method  demonstrated  in  previous  studies 
(Mlakar  1984)  (Mlakar  and  Toussi  September  1985,  October  1985)  (Bryant  and 
Mlakar  1987)  .  Two  characteristics  of  this  procedure  which  are  natural 
requirements  of  such  a  method  are:  (1)  the  procedure  rationally  addresses  the 
uncertainties  inherent  in  the  degradation  of  mass  concrete  due  to  freezing  and 
thawing,  and  (2)  the  procedure  is  mathematically  straightforward  for  implemen¬ 
tation  by  CE  offices. 

3.  Herein,  this  procedure  is  demonstrated  by  hindcast  application  to 
structural  features  at  CE  Civil  Works  structures.  Three  criteria  were  con¬ 
sidered  in  reviewing  potential  candidates  for  application  of  the  procedure. 
First,  an  appreciable  degree  of  measurable  damage  due  to  freezing  and  thawing 
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of  the  feature  was  necessary  to  make  a  quantifiable  comparison  with  the  hind- 
cast  prediction.  Next,  the  availability  of  data  required  for  application  of 
the  procedure,  e.g.  temperature  and  concrete  properties,  was  necessary  for 
rational  application.  Finally,  the  selected  candidate  should  be  representa¬ 
tive  of  typical  CE  projects.  Six  projects  were  initially  considered  as  poten¬ 
tial  case  study  candidates  -  Dashields  Locks  and  Dam,  Emsworth  Locks  and  Dam, 
Lower  Monumental  Lock  Wall,  Lockport  Lock,  and  Locks  7  and  8  on  the  Mononga- 
hela.  Two  candidates  which  best  met  these  criteria  were  selected  as  case 
studies  for  application  of  the  procedure.  These  two  candidates  were  the  mid¬ 
dle  wall  and  land  wall  at  Dashields  Lock. 

4.  The  procedure  for  predicting  service  life  of  mass  concrete  previ¬ 
ously  developed  is  reviewed  in  Part  II  of  this  report.  The  procedure  is  then 
applied  to  two  case  studies  to  demonstrate  the  procedure  and  determine  the 
quantitative  and  qualitative  results  for  real  structures  exhibiting  damage  due 
to  freezing  and  thawing.  In  Part  III,  the  procedure  is  applied  to  the  verti¬ 
cal  surface  of  the  middle  lock  wall  at  Dashields  Lock.  Application  of  the 
procedure  to  the  land  wall  at  Dashields  Lock  is  described  in  Part  IV.  The 
sensitivity  of  the  probabilistic  results  to  material  and  environmental  factors 
is  discussed  in  Part  V.  In  Part  VI,  current  and  proposed  testing  requirements 
relevant  to  the  prediction  of  concrete  service  life  are  discussed.  Finally, 
our  conclusions  and  recommendations  related  to  this  work  are  set  forth  in 
Part  VII. 
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PART  II:  PROBABILISTIC  PROCEDURE  FOR  PREDICTING  SERVICE  LIFE 


Prediction  Procedure 

5.  A  probabilistic  procedure  for  predicting  service  life  of  mass  con¬ 
crete  subject  to  damage  due  to  freezing  and  thawing  was  previously  presented 
(Bryant  and  Mlakar  1988) .  The  mechanism  of  freezing  and  thawing  used  in  this 
procedure  is  based  on  the  thesis  that  damage  in  a  unit  volume  occurs  when  the 
temperature  is  below  a  critical  value  (8  <  ©cr)  and  the  degree  of  saturation 
(amount  of  freezable  water)  exceeds  a  critical  value  (S  >  Scr) . 1  In  a  deter¬ 
ministic  sense,  this  event  is  depicted  in  Figure  1,  which  indicates  hypotheti¬ 
cal  time-histories  of  temperature  and  saturation  in  a  unit  volume  of  concrete. 
As  seen  in  this  figure,  the  temperature  and  degree  of  saturation  vary  in 
response  to  external  environmental  conditions  as  determined  by  the  material 
properties.  When  either  event  occurs  singly  there  is  no  damage  but  rather  it 
is  their  simultaneous  occurrence  that  causes  damage.  In  this  context,  sound¬ 
ness  of  the  concrete  is  considered  a  binary  state  process,  i.e.  either  the 
unit  volume  is  damaged  or  it  is  not. 

6.  The  previous  description  of  damage  due  to  freezing  and  thawing  is 
rational  but  deterministic.  In  fact,  the  temperature  and  degree  of  saturation 
in  a  structure  vary  randomly  in  time  due  to  randomness  of  environmental  condi¬ 
tions.  As  a  consequence  of  the  random  nature  of  these  conditions,  the  damage 
due  to  freezing  and  thawing  and  thus  service  life  is  also  a  random  variable. 

In  addition,  the  relevant  material  properties  and  critical  values  of  tempera¬ 
ture  and  saturation  introduce  additional  uncertainty  in  the  prediction  or 
service  life.  Under  these  conditions,  the  hazard  function  or  instantaneous 
rate  of  failure  is: 

h(t)  =  P ( 8  <  8cr  \  S  >  Scr}  (1) 

that  is,  the  joint  probability  of  temperature  less  than  critical  and  satura¬ 
tion  greater  than  critical.  Although  the  hazard  function  generally  is  a  func¬ 
tion  of  time,  on  the  scale  of  damage  due  to  freezing  and  thawing  it  is 

1  It  should  be  noted  that  damage  normally  will  occur  only  with  non-air- 
entrained  concrete. 
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reasonable  to  discuss  annual  failure  rates,  A  ,  which  are  unaffected  by  sea¬ 
sonal  variations.  For  simplicity  herein,  we  assume  that  the  factors  which 
determine  criticality  of  temperature  are  independent  of  those  defining  the 
criticality  of  saturation  such  that: 

X  -  P(8  <  0cr)  •  P { S  >  Scr)  (2) 

7.  The  previous  study  demonstrated  that  the  constant  hazard  function  X 
implies  an  exponential  distribution  of  service  life.  Further,  the  expected 
service  life  is  the  reciprocal  of  the  failure  rate,  i.e  : 

T  -  1/A  (3) 

It  is  obvious  from  Equations  2  and  3  above  that  an  increase  in  the  probability 
of  temperature  less  than  critical  or  saturation  greater  than  critical 
increases  the  failure  rate  and  inversely  decreases  the  service  life.  The 
behavior  implied  by  the  above  equations  agrees  with  the  observed  deterioration 
history  at  typical  CE  concrete  structures,  i.e.  some  structures  deteriorate 
rapidly  and  need  early  repair  while  others  experience  very  slow  deterioration. 
These  extremes  are  indicated  as  A3  and  A3  ,  respectively,  in  Figure  2.  The 
primary  interest  of  the  Repair,  Evaluation,  Maintenance,  and  Rehabilitation 
Research  Program  is  with  intermediate  cases  with  gradual  but  significant  dam¬ 
age,  such  as  that  indicated  by  X2  in  Figure  2. 

8.  In  the  discussion  above,  it  is  apparent  that  the  probabilistic 
procedure  includes  four  basic  steps  prio1'  to  determination  of  the  required 
probabilities  of  temperature  and  saturation  criticality  for  the  unit  volume. 
These  steps  are: 

a.  Determine  critical  temperature,  8cr  . 

b.  Determine  critical  degree  of  saturation,  Scr  . 

c.  Determine  distribution  of  temperature,  9  . 

d.  Determine  distribution  of  degree  of  saturation,  S  . 

9.  The  critical  temperature  is  a  function  of  several  factors  which 
affect  the  freezability  of  water  in  the  concrete.  The  impurities  in  the 
water  reduce  the  freezing  point  below  normal  and,  as  any  part  of  the  water  in 
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a  pore  begins  to  freeze,  the  remaining  part  has  an  even  lower  freezing  point.* 
At  this  time  these  effects  have  not  been  quantified  and,  for  the  purposes  of 
this  study,  a  constant  deterministic  value  of  0cr  will  be  assumed. 

10.  The  critical  degree  of  saturation  is  determined  by  a  number  of 
factors  not  known  with  certainty.  Critical  pore  size,  moisture  migration  due 
to  freezing,  and  frost  generated  pressures  are  known  to  affect  the  critical 
degree  of  saturation.  Currently,  no  widely  accepted  quantitative  relationship 
between  these  variables  is  available.  The  previous  study  (Bryant  and  Mlakar 
1988)  suggested  that  a  degree  of  saturation  less  than  91  percent  is  not  criti¬ 
cal,  based  on  the  9  percent  expansion  of  water  when  frozen. 

11.  Temperature  in  a  unit  volume  can  be  estimated  from  the  environmen¬ 
tal  conditions  and  thermal  material  properties.  For  some  simpler  cases, 
closed-form  solutions  of  the  deterministic  heat  transfer  equations  exist 
(Carslaw  and  Jaeger  1959).  More  generally,  finite-element  or  finite- 
difference  computer  codes  provide  a  solution  (e.g.  ADINA  Engineering  1987). 

For  the  probabilistic  analysis,  we  need  to  determine  the  distribution  of 
temperature  for  the  unit  volume,  as  illustrated  in  Figure  3.  In  determination 
of  this  distribution,  the  uncertainty  or  randomness  of  the  functionally  inde¬ 
pendent  variables  in  the  deterministic  analysis  must  be  considered.  In  gen¬ 
eral,  these  random  variables  include  environmental  temperature  0g  ,  thermal 
conductivity  k  ,  specific  heat  Cp  ,  and  material  density  p  . 

12.  Calculation  of  the  degree  of  saturation  experienced  by  a  unit 
volume  of  concrete  presents  a  more  challenging  problem,  even  in  the  determin¬ 
istic  sense.  The  previous  study  indicated  that  moisture  migration  in  concrete 
is  generally  considered  to  include  two  mechanisms  -  seepage  and  capillary 
action.**  Seepage  through  permeable  materials  is  reasonably  well  understood 
and  depends  primarily  on  the  coefficient  of  permeability  of  the  material  and 
the  external  heads  of  the  fluid.  On  the  other  hand,  capillary  action  is 
understood  only  conceptually,  with  little  actual  measurements  or  predictions 


*  Assuming  the  water  is  not  a  saturated  solution  before  freezing  begins  -  as 
it  will  be  if  nothing  is  dissolved  in  it  but  calcium  hydroxide. 

**  In  fact;  for  typical  concrete  having  any  value  of  water-cement  ratio  below 
0.7,  there  will  be  no  capillary  continuity  in  concrete  after  1  vear;  at 
lower  valves  continuity  will  be  lost  sooner,  by  14  days  for  w/c  -  0.5, 
hence  moisture  movement  will  be  by  surface  diffusion  in  intact  concrete. 
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available.  For  the  purposes  of  these  case  studies,  moisture  migration,  and 
thus  degree  of  saturation,  is  computed  only  by  considering  the  seepage  mecha¬ 
nism.  For  the  probabilistic  analysis,  we  determine  the  distribution  of  degree 
of  saturation  for  the  unit  volume,  as  illustrated  in  Figure  4.  In  determina¬ 
tion  of  this  distribution,  the  randomness  of  the  functionally  independent 
variables  in  the  deterministic  analysis,  permeability  and  external  head,  are 
considered. 


Environmental  Data 


13.  The  environmental  data  required  for  prediction  of  damage  due  to 
freezing  and  thawing  by  this  procedure  includes  air  temperatures,  water  tem¬ 
peratures,  ground  temperatures,  and  water  levels.  Each  of  these  input  random 
variables  is  described  by  a  probability  distribution  function.  The  data  for 
each  variable  was  determined  from  an  appropriate  source  for  the  case  study 
site  and  represented  by  its  statistical  parameters. 

14.  For  simplification  of  the  input  to  the  thermal  analyses,  the 
annual  variations  of  air  and  water  temperature  were  represented  by  a  sinusoi¬ 
dal  variation  of  temperature  with  an  annual  mean  and  amplitude.  These  varia¬ 
tions  are  defined  as: 

0a  —  A  +  B  sin(ut)  +  C  sin(wt)  (4) 

-  A  +  D  sin(wt  +  <j>) 

0W  -  E  +  F  sin(cot)  +  G  sin(wt)  (5) 

-  E  +  H  sin(u>t  +  4>) 

where : 

0a  -  Air  temperature. 

0W  -  Water  temperature . 
w  -  Frequency  -  2p/1  . 

4>  -  Phase  shift. 

and  A  through  H  are  constants  determined  by  linear  regression  of  a  sinusoi¬ 
dal  fit  to  the  data,  as  described  in  Appendix  A.  An  example  of  such  data  is 
in  Table  1 . 
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PART  III:  CASE  STUDY  OF  MIDDLE  WALL  AT  DASHIELDS  LOCK 


15.  The  first  case  study  focused  on  the  middle  wall  (between  the  two 
locks)  at  Dashields  Lock  on  the  Ohio  River  near  Pittsburgh.  The  layout  of 
Dashields  Lock  is  shown  in  Figure  5,  with  a  section  through  the  middle  wall 
shown  in  Figure  6.  Key  features  of  the  middle  wall  include  an  access  gallery 
in  the  top  section  and  the  filling  and  emptying  gallery  in  the  lower  section. 
The  wall  is  exposed  to  water  in  the  lock  chambers  on  both  sides. 

Temperature  Criticality  -  Middle  Wall 

16.  The  air  temperature  data  were  obtained  from  Local  Climatological 
Data  provided  by  the  National  Oceanic  and  Atmospheric  Administration  (1986) . 
For  Dashields  Lock,  the  air  temperatures  compiled  for  Pittsburgh,  PA,  were 
used  due  to  its  close  proximity  to  the  site.  Average  monthly  air  temperatures 
for  Pittsburgh,  shown  in  Figure  7,  were  the  basis  for  the  hindcast 
prediction . 

17.  Water  cemperature  data  at  the  site  were  obtained  from  the  Ohio 
River  Sanitary  Commission.  Water  temperature  data  covering  4  years  of  daily 
measurement  were  downloaded  from  the  STORET  data  base  by  US  Army  Engineer 
Waterways  Experiment  Station  (WES)  personnel.  These  data  were  synthesized 
into  an  annual  temperature  variation,  shown  in  Figure  8,  by  averaging  water 
temperatures  of  each  calendar  day  over  the  4  years. 

18.  Data  were  not  available  for  annual  variation  of  ground  temperatures 
at  the  site.  It  was  assumed  that  the  ground  temperature  at  the  base  of  the 
lock  could  be  reasonably  represented  as  constant,  and  a  value  of  50°  F  was 
selected. 

19.  The  temperature  boundary  conditions  for  thermal  analysis  of  the 
middle  wall  are  indicated  in  Figure  9.  Above  the  mean  waterline,  the  surface 
temperature  was  assumed  to  vary  in  sinusoidal  fashion  as  shown  in  Figure  10. 
The  temperature  at  the  wall  surface  below  the  waterline  was  constrained  to 
vary  as  shown  in  Figure  11.  At  the  base  of  the  structure,  the  ground  temper¬ 
ature  was  assumed  to  remain  constant  at  50°  F.  The  upper  gallery  was  assumed 
open  to  the  air  and  thus  exposed  to  the  air  temperatures  in  Figure  10.  The 
lower  gallery,  continuously  filled  with  water,  was  therefore  constrained  to 
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water  temperatures.  The  sinusoidal  fits  to  the  air  and  water  temperatures 
according  to  Equations  4  and  5  produced  the  constants  summarized  in  Table  1 . 
The  comparisons  in  Figures  10  and  11  confirm  that  this  fit  is  a  good  repre¬ 
sentation  of  the  actual  data. 

20.  Calculation  of  temperatures  for  a  unit  volume  of  concrete  was 
accomplished  by  using  two  analysis  techniques,  i.e..  one -dimensional  (1-D) 
and  two-dimensional  (2-D)  models.  The  1-D  analysis  uses  an  analytical  solu¬ 
tion  of  a  slab  with  periodic  temperature  variation  on  each  surface  after 
Carslaw  and  Jaeger  (1959)  as  indicated  in  Figure  12.  For  the  periodic  varia¬ 
tion  of  external  temperature  0e  ,  the  attenuation  of  temperature  with  dis¬ 
tance  from  the  surface  is  predicted  by  the  function  X(x)  shown  in  the  fig¬ 
ure.  As  evident  from  this  expression,  temperature  is  a  function  of  external 

O 

temperature  0e  ,  distance  from  the  surface  x  ,  and  thermal  diffusivity  h'1  . 
This  attenuation  of  temperature  with  distance  is  indicated  (for  a  unit  0e)  in 
Figure  13  for  three  typical  values  of  thermal  diffusivity. 

21.  Two  1-D  thermal  analyses  of  the  wall  were  conducted  for  locations 
2  ft  above  and  below  the  waterline.  The  attenuation  of  temperature  with 
depth,  previously  indicated  in  Figure  13,  also  indicates  the  sensitivity  to 
thermal  diffusivity.  The  annual  temperature  cycles  above  the  waterline  at 
depth  intervals  of  2  ft*  are  shown  in  Figure  14.  This  figure  presents  the 
same  information  as  Figure  13  in  a  temporal  format.  An  instructive  look  at 
calculation  of  the  temperature  distribution  for  a  depth  of  6  in.  is  presented 
in  Appendix  B.  These  sample  calculations  demonstrate  the  procedure  for  cal¬ 
culation  of  the  deterministic  temperatures  and  Taylor- series  expansion  for 
temperature  distribution.  A  most  important  fact  exposed  by  the  calculation  of 
the  temperature  distribution  is  the  relative  contribution  of  the  various  fac¬ 
tors  to  the  variance  of  temperature.  The  calculations  indicate  that,  in  this 
case,  over  99  percent  of  the  overall  variance  is  due  to  the  annual  temperature 
amplitude,  0a  .  This  important  factor  was  verified  for  all  depths  at  the  two 
locations  above  and  below  the  waterline  as  shown  in  Appendix  B.  In  all  cases, 
the  temperature  amplitude  accounted  for  over  97  percent  of  the  overall  vari¬ 
ance.  Thus,  the  uncertainties  in  other  factors  can  be  rationally  ignored  in 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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further  thermal  analyses.  This  conclusion  permits  estimation  of  the  overall 
temperature  distribution  with  only  a  single  deterministic  analysis  at  the  mean 
values,  an  extremely  helpful  factor  in  the  2-D  analyses. 

22.  The  2-D  finite -element  model  of  the  WES  Concrete  Technology  Divi¬ 
sion  (CTD)  for  the  thermal  analysis  using  the  ABACUS  program  is  illustrated  in 
Figure  15.  The  symmetry  of  geometry  and  boundary  conditions  permitted  analy¬ 
sis  of  only  one-half  of  the  wall,  with  a  constraint  of  no  heat  flow  across  the 
plane  of  symmetry.  The  mesh  is  more  highly  discretized  in  the  region  of  most 
interest  and  highest  temperature  variations,  that  is,  along  the  wall  surface 
near  the  waterline.  A  transient  thermal  analysis  was  conducted  for  a  period 
of  3  years  to  obtain  an  annual  cycle  of  temperatures  at  all  nodal  points  in 
the  wall.  The  1-D  analytical  solution  indicated  that  the  effect  of  initial 
conditions  would  disappear  after  less  than  1  year.  Thus,  the  third  year 
annual  cycle  was  used  to  characterize  the  temperature  distribution  at  any 
point  in  the  wall. 

23.  At  each  location  in  both  the  1-D  and  2-D  analyses,  the  mean  and 
annual  amplitude  of  temperature  was  obtained  by  fitting  a  sinusoidal  curve  to 
the  computed  annual  temperature  cycle.  This  curve  fit  was  accomplished  in 
the  same  manner  as  described  in  Appendix  A  for  the  environmental  temperatures. 
An  example  of  the  curve  fit  is  illustrated  in  Figure  16,  which  is  shown  for  a 
typical  node.  The  annual  cycle  of  temperatures  and  the  sinusoidal  fit  are 
shown  to  agree  very  well.  A  comparison  of  mean  temperatures  and  annual  ampli¬ 
tudes  from  the  1-D  and  2-D  analyses  are  compared  in  Figures  17  and  18  for  the 
locations  2  ft  above  and  below  the  waterline,  respectively.  As  seen  in  these 
comparisons,  the  1-D  analyses  predict  reasonable  values  for  both  mean  and 
amplitude,  with  the  differences  attributable  to  the  actual  2-D  nature  of  the 
problem.  For  example,  above  the  waterline  (Figure  17)  the  effect  of  the  tem¬ 
peratures  below  the  waterline  is  evidenced  in  the  2-D  results  but  not  in  the 
1-D  results.  The  proximity  of  the  warmer  temperatures  increases  the  mean 
temperature  away  from  the  air  temperature  boundary.  The  gallery  reduces  the 
annual  temperature  amplitude  computed  in  the  2-D  analysis.  In  effect,  the 
gallery  has  the  same  effect  as  if  the  wall  were  actually  less  than  24  ft 
thick.  Similar  effects  are  seen  in  Figure  18  for  the  location  below  the 
waterline . 


24.  The  annual  probability  of  freezing  at  these  two  locations,  accord¬ 
ing  to  the  1-D  and  2-D  analyses,  is  compared  in  Figure  19.  Above  the  water¬ 
line,  the  probabilities  range  from  about  20  percent  at  the  surface  to  zero  at 
2.5  ft.  Below  the  waterline,  the  temperatures  never  fall  below  32°  F.  These 
comparisons  instill  confidence  in  the  2-D  results  and  indicate  that  in  some 
cases,  1-D  analyses  provide  reasonable  results. 

25.  The  results  of  the  2-D  thermal  analysis  are  presented  in  detail 
for  the  area  of  the  wall  indicated  in  Figure  20.  The  computed  mean  tempera¬ 
ture  levels  are  contoured  in  Figure  21.  As  expected,  the  mean  temperatures 
near  the  surface  reflect  the  mean  air  and  water  temperatures  and  transition 
between  these  extremes.  The  annual  amplitudes  of  temperature  are  depicted  for 
this  area  in  Figure  22.  Near  the  surface,  the  amplitudes  are  near  those  of 
the  air  and  water  temperatures,  which  are  almost  the  same.  In  addition  to  the 
expected  reduction  of  amplitude  with  distance  from  the  surface,  the  effect  of 
the  upper  and  lower  galleries  is  evident  in  this  plot.  As  the  surfaces  of  the 
gallery  walls  are  exposed  to  the  same  temperature  variations  as  the  outer 
surfaces,  the  temperature  amplitudes  near  these  surfaces  tend  to  have  higher 
values.  Probability  of  freezing  in  this  area  of  the  wall  is  illustrated  in 
Figure  23.  As  should  be  expected,  the  probability  of  freezing  is  zero  below 
the  waterline  and  for  some  interior  portions  above  the  waterline.  The  proba¬ 
bility  of  freezing  above  the  waterline  is  primarily  driven  by  the  reduction  of 
temperature  amplitude  with  increasing  depth. 

Saturation  Criticality  -  Middle  Wall 

26.  As  noted  in  Part  II,  determination  of  degree  of  saturation  in  con¬ 
crete  is  a  complex  problem  considered  to  embody  two  methods  of  moisture 
migration  -  seepage  flow  and  capillary  action.  Although  there  have  been 
attempts  in  the  past  to  estimate  height  of  capillary  rise  in  various  soils, 
these  estimates  are  extremely  approximate  and  do  not  extend  to  other  materials 
such  as  concrete.  For  the  purposes  of  these  case  studies,  moisture  migration 
was  considered  to  occur  by  seepage  alone.  Further,  although  2-D  finite- 
element  and  finite-difference  solutions  for  seepage  exist,  a  simpler  solution 
to  this  essentially  steady-state  problem  was  applied.  This  approach,  depicted 
in  Figure  24,  reflects  a  steady-state  location  of  the  phreatic  surface  in  a 
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surface  in  a  concrete  wall  for  a  given  set  of  boundary  conditions  (external 
head).  Essentially,  in  this  case,  the  phreatic  surface  is  along  a  straight 
line  between  the  water  levels  on  each  side  of  the  wall.  Attendant  to  this 
assumption  is  the  recognition  that  the  material  is  fully  saturated  along  and 
below  the  phreatic  surface  and  unsaturated  above  it.  The  seepage  model  gives 
no  information  on  the  expected  variation  of  degree  of  saturation  above  the 
phreatic  surface.  Thus,  the  location  of  91  percent  saturation  cannot  be 
determined  from  the  seepage  model.  Therefore,  for  these  case  studies,  the 
degree  of  saturation  is  considered  critical  below  the  phreatic  surface  and  not 
critical  above. 

27.  The  water  levels  in  a  lock  were  characterized  by  a  mixed  distribu¬ 
tion  shown  in  Figure  25.  This  distribution  is  based  on  normal  distributions 
with  means  at  high  and  low  pool  elevations  of  690  and  684  ft,  respectively , 
with  a  standard  deviation  about  each  mean  of  6  in.  The  combined  distribution 
assumes  equal  probability  of  high  or  low  pool  and  is  adjusted  to  maintain  a 
maximum  probability  of  unity.  As  the  water  levels  on  each  side  of  the  wall 
have  the  same  distribution,  the  location  of  the  phreatic  surface,  and  thus  the 
location  of  critical  saturation,  has  the  same  distribution  as  the  water  level. 
Attendant  to  this  observation,  the  probability  of  critical  saturation  is  as 
shown  in  Figure  26.  This  plot  reflects  the  cumulative  distribution  function 
(CDF)  of  Figure  25.  It  should  be  noted  from  the  CDF  that  the  probability  of 
critical  saturation  is  a  constant  50  percent  from  about  elevation  685  to 
689  ft. 


Material  Service  Life  -  Middle  Wall 


28.  As  a  preliminary  step  in  calculation  of  the  service  life,  the 
annual  hazard  function  is  computed  as  described  by  Equation  2.  This  hazard 
function  is  depicted  for  the  central  section  of  the  middle  wall  in  Figure  27. 
As  defined  by  Equation  2,  this  figure  is  the  product  of  Figures  23  and  26. 

As  should  be  expected,  the  highest  probability  of  damage  (the  hazard  function) 
is  near  the  wall  surface  just  above  the  waterline.  This  is  the  area  of  higher 
probabilities  of  critical  saturation  and  critical  temperatures.  Naturally, 
the  area  below  the  waterline  and  near  the  interior  of  the  wall  has  essentially 
no  chance  of  damage  since  the  temperatures  have  little  chance  of  being  below 
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critical  in  these  areas.  As  seen  in  Figure  26,  the  probability  of  saturation 
is  nil  more  than  2  ft  above  the  mean  high  pool,  thus  the  hazard  function  is 
likewise  nil  in  this  area. 

29.  The  material  service  life  for  degradation  due  to  freezing  and  thaw¬ 
ing  is  estimated  for  the  middle  wall  according  to  Equation  3.  This  function, 
the  inverse  of  the  hazard  function,  is  plotted  in  Figure  28.  As  seen  in  this 
plot,  the  service  life  is  predicted  to  be  only  15  years  at  a  depth  of  9  in. 

At  2  ft,  the  service  life  is  predicted  to  be  about  40  years,  with  a  service 
life  of  about  100  years  at  3  ft. 

30.  This  service  life  prediction  (hindcast,  actually)  is  encouragingly 
similar  to  the  observed  deterioration  on  the  middle  wall  at  this  site.  Mea¬ 
surements  of  damage  due  to  freezing  and  thawing  in  the  middle  wall  made  by  the 
WES  (1986)  in  1980  indicated  depths  of  deterioration  of  0  to  2  ft.  Given  the 
fact  that  the  lock  was  constructed  in  1927-1929,  the  measurements  of  damage 
after  50  years  are  very  close  to  those  hindcast  by  this  procedure,  although 
damage  rates  on  this  structure  were  affected  by  at  least  two  additional  fac¬ 
tors.  One,  some  repair  work  (shotcrete  overlays)  would  tend  to  reduce  rate  of 
damage.  Two,  abrasion  by  barges  would  tend  to  enhance  the  rate  of  damage. 
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PART  IV:  CASE  STUDY  OF  LAND  WALL  AT  DASHIELDS  LOCK 


31.  The  second  case  study  was  the  land  wall  at  Dashields,  previously 
shown  in  Figure  5.  A  section  through  the  land  wall,  shown  in  Figure  29, 
serves  to  indicate  the  similarities  and  differences  of  the  land  wall  and  the 
middle  wall.  The  wall  is  of  similar  construction  with  a  maximum  thickness  of 
24  ft  and  having  similar  upper  and  lower  galleries.  The  wall  is  not  symmet¬ 
ric,  stepping  down  to  a  minimum  thickness  of  12  ft  in  the  upper  section.  The 
other  primary  difference  between  the  two  walls  lies  in  the  boundary  condi¬ 
tions.  Obviously,  the  river  side  of  the  land  wall  is  exposed  to  similar  tem¬ 
perature  and  moisture  conditions  as  the  middle  wall.  The  land  side  of  the 
wall,  being  in  full  contact  with  the  ground,  is  exposed  to  different  tempera¬ 
ture  and  moisture  conditions. 


Temperature  Criticality  -  Land  Wa1  T 


32.  The  temperature  boundary  conditions  for  thermal  analysis  of  the 
land  wall  are  indicated  in  Figure  30.  Along  the  river  side  of  the  wall  and 
the  upper  and  lower  surfaces,  the  temperature  boundary  conditions  shown  are 
the  same  as  for  the  middle  wall.  Temperatures  in  the  upper  and  lower  gal¬ 
leries  were  likewise  constrained  to  the  ambient  air  and  water  temperatures,  as 
in  the  analysis  of  the  middle  wall.  The  temperatures  along  the  land  side  of 
the  wall  obviously  must  vary  between  the  annual  air  temperature  cycle  at  the 
surface  and  the  constant  ground  temperature  at  the  base.  In  order  to  estimate 
this  variation,  a  1-D  analysis  of  the  ground  as  a  semi - infinite  medium  exposed 
to  cyclic  temperature  variation  at  the  surface  was  performed  as  described  by 
Carslaw  and  Jaeger  (1959) .  This  1-D  analysis  produced  the  results  shown  in 
Figure  31.  These  results  indicate  that  the  cyclic  temperature  effect  has  an 
effect  only  for  about  the  top  5  ft  of  the  ground.  Below  this  depth,  the 
nature  of  the  cyclic  air  temperature  has  essentially  no  effect.  For  compara¬ 
tive  purposes,  the  results  are  compared  with  those  for  a  typical  concrete. 

33.  In  the  finite-element  thermal  analysis  by  the  CTD,  a  portion  of 
the  soil  was  modeled  in  an  attempt  to  better  define  the  appropriate  thermal 
conditions  along  the  land  face  of  the  wall.  The  finite-element  model,  shown 
in  Figure  32,  included  a  representation  of  the  soil  out  to  a  distance  of 
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30  ft  from  the  wall  at  the  ground  surface  (18  ft  at  the  base).  As  seen  in 
Figure  32,  an  adiabatic  boundary  (no  heat  flow)  was  specified  vertically  to 
an  elevation  of  680  ft  with  constant  ground  temperatures  (50°  F)  specified 
for  the  remainder  of  the  depth  Ambient  air  temperatures  were  specified 
along  the  ground  surface .  This  soil  portion  of  the  model  was  felt  to  better 
represent  actual  temperature  gradients  along  the  land  face  of  the  wall.  The 
results  near  the  river  face  of  the  wall,  where  temperature  and  saturation 
criticality  were  most  important,  are  not  likely  to  be  much  affected  by  this 
boundary  condition. 

34.  The  thermal  analysis  was  conducted  by  the  WES  CTD  using  the  ABACUS 
computer  code  for  the  land  wall  given  the  above  boundary  conditions.  In  this 
case,  the  entire  model  was  analyzed,  since  there  was  no  symmetry  as  in  the 
middle  wall.  As  for  the  analysis  of  the  middle  wall,  a  transient  analysis 
for  3  years  was  conducted,  and  the  nodal  temperature  results  of  the  third 
year  were  used  to  characterize  the  sinusoidal  distribution  of  internal  temper¬ 
atures.  This  analysis  produced  the  contours  of  mean  temperature,  shown  in 
Figure  33,  and  the  annual  variation  of  temperature,  shown  in  Figure  34,  for 
the  central  portion  of  the  lock  wall.  As  was  seen  in  the  middle  wall,  the 
mean  temperatures  and  annual  amplitudes  track  the  ambient  conditions  near 
each  surface  and  transition  rather  smoothly  in  this  section.  The  probability 
of  temperature  being  critical,  for  simplicity  of  agreement  taken  to  be  below 
32°  F,  varies  from  zero  to  about  0.18  in  this  region,  and  is  plotted  in  Figure 
35.  The  probability  of  freezing  exists  only  above  the  waterline,  and  gene¬ 
rally  decreases  with  distance  from  a  surface  exposed  to  the  air. 

Saturation  Criticality  -  Land  Wall 

35.  The  saturation  model  for  the  land  wall  is  quite  similar  to  that 
for  the  middle  wall,  recognizing  the  difference  in  water  levels  on  both  sides 
of  the  lock  wall.  As  for  the  middle  wall,  the  water  levels  in  the  lock  were 
characterized  by  the  mixed  distribution  shown  in  Figure  25.  The  actual  dis¬ 
tribution  of  water  table  elevation  on  the  land  side  of  this  wall  is  not 
known.  However,  piezometer  data  was  available  for  three  locations  along  the 
land  wall  for  a  period  of  about  4  years.  This  transient  data  for  the  piezo¬ 
meter  near  the  middle  of  the  lock  is  plotted  in  Figure  36  for  this  period. 
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For  the  purposes  of  the  case  study,  it  was  reasonable  to  assume  that  the  land 
side  water  table  fluctuated  much  more  slowly  than  the  lock  pool  elevation  and 
that  it  could  be  reasonably  represented  by  an  average  water  table  elevation, 
also  indicated  as  about  690  ft  in  Figure  36.  As  the  phreatic  surface  was 
assumed  to  lie  along  a  line  between  the  water  levels  on  each  side  of  the  wall 
(as  for  the  middle  wall),  the  resulting  distribution  of  this  surface,  and  thus 
the  probability  of  critical  saturation,  is  as  depicted  in  Figure  37.  Essen¬ 
tially,  the  distribution  at  the  river  side  of  the  wall  (shown  in  Figure  25)  is 
mapped  through  the  depth  of  the  wall,  decreasing  to  the  point  distribution  at 
the  land  face  of  the  wall.  Thus,  the  distribution  at  any  distance  from  the 
surface  is  identical,  but  acts  over  a  smaller  range  of  elevation.  The  proba¬ 
bility  of  critical  saturation  was  computed  at  every  point  in  the  wall  corre¬ 
sponding  to  the  nodal  points  in  the  thermal  analysis,  and  this  probability  is 
plotted  in  Figure  38.  This  contour  plot  reflects  the  mapping  of  the  pool 
elevation  distribution  at  the  river  side  to  the  constant,  mean  water  table 
elevation  at  the  land  side. 

Material  Service  Life  -  Land  Wall 

36.  The  annual  hazard  function,  as  defined  by  Equation  2,  is  depicted 
for  the  central  section  of  the  land  wa] 1  in  Figure  39.  Once  again,  the  high¬ 
est  probability  of  damage  is  near  the  wall  surface  just  above  the  water-line. 
This  is  the  area  of  higher  probabilities  of  critical  saturation  and  critical 
temperatures.  The  area  below  the  waterline  and  near  the  interior  of  the  wall 
has  essentially  no  chance  of  damage  since  the  temperatures  have  little  chance 
of  being  below  critical  in  these  areas.  As  was  seen  in  Figure  38,  the  zero 
probability  of  saturation  more  than  2  ft  above  the  mean  high  pool  produces  a 
zero  probability  of  damage  in  this  area. 

37.  The  material  service  life  for  degradation  due  to  freezing  and  thaw¬ 
ing  for  the  land  wall  is  estimated  according  to  Equation  3,  and  is  plotted  in 
Figure  40.  As  seen  in  this  plot,  the  service  life  is  predicted  to  be  only  15 
years  at  a  depth  of  about  6  in.,  with  service  life  about  75  years  at  2  ft. 

38.  This  service  life  prediction  is  generally  similar  to  the  observed 
deterioration  on  the  land  wall  at  this  site.  Specifically,  measurements  of 
damage  due  to  freezing  and  thawing  in  the  land  wall  (WES  1986)  in  1980 
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indicated  depths  of  deterioration  generally  somewhat  less  than  along  the  mid¬ 
dle  wall.  The  constraining  effect  of  the  relatively  constant  ground  tempera¬ 
tures  and  the  thinner  upper  section  of  the  wall  reduced  the  area  subject  to 
critical  temperatures  in  the  land  wall.  As  a  result,  the  probability  of  dam¬ 
age  is  reduced,  and  the  expected  service  life  increased  in  the  land  wall. 
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PART  V:  SENSITIVITY  TO  KEY  MATERIAL  AND  ENVIRONMENTAL  VARIABLES 


39.  The  calculations  presented  in  Appendix  B  provided  an  instructive 
look  at  calculation  of  the  Taylor-series  expansion  for  temperature  distribu¬ 
tion.  An  important  fact  revealed  by  this  calculation  of  the  temperature  dis¬ 
tribution  was  the  relative  contribution  of  the  various  factors  to  the  variance 
of  temperature.  The  calculations  indicated  that,  for  a  depth  of  6  in.,  over 
99  percent  of  the  overall  variance  is  due  to  the  annual  temperature  amplitude. 
In  fact,  the  temperature  amplitude  accounted  for  over  97  percent  of  the  over¬ 
all  variance  for  all  cases  evaluated  in  the  1-D  analyses.  This  fact  indicated 
that  the  uncertainties  in  other  factors  could  be  rationally  ignored  in  the 
thermal  analyses.  This  conclusion  permitted  estimation  of  the  overall  temper¬ 
ature  distribution  with  a  single  deterministic  analysis  at  the  mean  values. 

In  terms  of  sensitivity,  the  temperature  distribution  is  sensitive  only  to  the 
annual  amplitude  of  temperature.  Uncertainty  in  other  factors  did  not  sig¬ 
nificantly  contribute  to  uncertainty  of  internal  concrete  temperatures. 


Saturation  Distribution 

40.  The  estimation  of  the  degree  of  saturation  in  the  concrete  was 
based  on  three  general  assumptions.  First,  the  water  level  in  the  lock  was 
assumed  to  follow  the  distribution  previously  indicated  in  Figure  2b.  This 
mixed  distribution  assumes  two  normal  distributions  centered  about  the  nominal 
high  and  low  pool  elevations,  which  are  the  nominal  high  and  low  water  eleva¬ 
tions  in  the  lock.  A  reasonable  standard  deviation  about  these  two  means  was 
assumed  to  represent  fluctuations  in  the  two  pool  elevations.  This  distribu¬ 
tion  reflects  the  assumption  that  the  water  in  the  lock  is  left  at  whatever 
level  exists  when  the  exit  gate  is  closed  until  the  next  lockage  occurs. 

Since  there  was  no  evidence  to  the  contrary,  it  was  assumed  that  there  was  an 
equal  probability  of  downstream  or  upstream  traffic,  and  thus  the  two  normal 
distributions  are  weighted  equally. 

41.  The  second  general  assumption  in  development  of  the  degree  of 
saturation  involved  the  method  of  moisture  migration  and  the  procedure  used  to 


determine  the  degree  of  saturation.  Moisture  migration  was  assumed  to  occur 
primarily  by  seepage.  For  the  problem  at  hand  with  small  differences  in 
external  heads,  the  phreatic  surface  lies  essentially  on  a  straight  line 
between  the  external  water  levels,  and  the  straight  line  surface  was  utilized 
in  the  analysis. 

42.  The  third  assumption  in  determination  of  saturation  criticality 
was  that  the  phreatic  surface  (100  percent  saturation)  defined  the  boundary 
of  critical  saturation.  The  seepage  analysis  for  unconfined  flow  provides 
only  this  surface  and  does  not  predict  where  a  lesser  degree  of  saturation, 
such  as  91  percent,  exists.  For  the  purposes  of  this  study,  this  effect  was 
felt  to  be  small. 

43.  For  the  assumptions  of  this  study,  the  primary  uncertainty  in 
degree  of  saturation  and  probability  of  critical  saturation  resides  totally 
in  the  assumed  distributions  of  external  head.  Other  uncertainties  not  con¬ 
sidered  in  this  study  might  be  as  important  as  external  head  and  are  worthy  of 
additional  study. 
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PART  VI:  TESTING  FOR  SERVICE  LIFE 


44.  The  procedure  presented  herein  for  prediction  of  concrece  service 
life  presumes  the  availability  of  information  regarding  environmental  condi¬ 
tions,  thermal  and  permeability  properties  of  concrete,  and  methods  to  compute 
the  response  of  the  structure  for  these  conditions.  In  addition  to  knowing 
expected  or  mean  values  of  the  input  variables,  the  uncertainty  associated 
with  each  is  also  needed.  The  current  availability  of  this  information  indi¬ 
cates  that  physical  testing  is  warranted  to  improve  the  knowledge  in  some 
areas . 

45.  Experimental  work  supporting  prediction  of  concrete  service  life 
obviously  includes  measurement  of  temperature  and  moisture  content  of  concrete 
by  both  laboratory  and  in-situ  testing.  Experimental  work  in  these  areas  can 
provide  valuable  information  in  two  different,  but  related  forms.  First, 
physical  tests  can  provide  first-hand  information  and  data  regarding  the 
actual  conditions  at  a  site  for  evaluation  of  the  condition  of  the  concrete. 
Second,  and  important  for  the  long-term  improvement  of  service  life  predic¬ 
tion,  these  measurements  can  provide  validation  for  analytical  methods  of 
temperature  and  moisture  prediction. 


Temperature 


46.  Measurements  of  temperature  for  prediction  of  concrete  service 
life  would  be  most  appropriately  obtained  by  in-situ  testing.  Temperature 
measurements  are  required  both  for  the  external  conditions,  such  as  air, 
water,  and  ground  temperatures,  and  internal  to  the  structure  itself.  These 
measurements  should  be  made  on  a  long-term  basis,  i.e.  over  a  period  of 
years,  to  evaluate  the  external  temperatures  and  the  thermal  response  of  the 
structure.  Current  procedures  and  measurement  devices  are  certainly  adequate 
to  provide  this  information.  The  major  difficulty  in  thermal  measurement  lies 
in  the  difficulty  in  delivery  or  placement  of  fucoe  devices  internal  to  exist¬ 
ing  structures.  Methods  to  place  thermocouples  in  a  structure  with  minimal 
cost  and  minimal  disturbance  of  the  structure  and  its  function  should  be  the 
primary  concern  at  this  time.  The  measurement  and  recording  devices  should 
also  provide  sufficient  reliability  to  provide  accurate  measurements  over  the 
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period  of  years  required  to  accumulate  the  necessary  data.  Analytical  studies 
of  the  thermal  response  of  the  structure  should  be  an  integral  part  of  the 
test  planning  to  insure  appropriate,  strategic  placement  of  gages. 


Saturation 


47.  Evaluation  of  degree  of  saturation  in  concrete  should  be  considered 
in  both  in- situ  and  laboratory  tests.  Core  samples  removed  from  a  subject 
structure  can  be  sealed  to  minimize  moisture  migration  and  tested  in  a 
laboratory.  A  current  American  Society  for  Testing  and  Materials  (ASTM)  pro¬ 
cedure,  ASTM  C  642-82  (ASTM  1989)  could  provide  a  basis  for  laboratory  deter¬ 
mination  of  degree  of  saturation  in  concrete.  This  procedure  currently 
addresses  measurements  for  specific  gravity,  absorption,  and  voids  in  hardened 
concrete.  Two  slight  modifications  of  the  procedure  would  allow  estimation  of 
in-situ  degree  of  saturation.  First,  the  specimen  should  be  sealed  with  a 
relatively  impermeable  coating  immediately  after  removal  from  the  structure  to 
prevent  drying  or  wetting  of  the  concrete.  After  this  coating  is  removed  in 
the  lab,  the  mass  of  the  specimen  should  be  recorded.  The  remainder  of  the 
procedure,  including  complete  drying  and  saturation  of  the  specimen,  would 
remain  as  specified  in  the  standard.  The  in-situ  degree  of  saturation  could 
then  be  computed  from  the  original  mass,  dry  mass,  and  fully  saturated  mass 
as : 


s  -  <w0  -  Vd)/(WS  -  Wd) 


(6) 


where : 

S  -  Degree  of  saturation,  percent. 

Wq  -  Original  specimen  mass. 

Wd  -  Dry  specimen  mass. 

Ws  -  Saturated  specimen  mass. 

This  modified  procedure  should  provide  a  reasonable  estimation  of  the  actual 
degree  of  saturation  in  the  structure  at  the  location  of  the  specimen. 

48.  In-situ  measurements  for  degree  of  saturation  are  expected  to  pro¬ 
vide  considerably  more  challenges  than  those  for  temperature.  Measurements  of 
external  factors,  such  as  water  levels  and  rainfall  amounts,  should  be 
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attainable  using  currently  available  technology.  The  major  consideration  in 
these  measurements  would  be  the  long-term  reliability  of  the  measuring  and 
recoiding  devices.  A  first  cut  at  measurement  of  saturation  in  a  concrete 
structure  would  probably  reflect  the  current  technology  in  measurement  of 
moisture  in  other  media  such  as  soil.  That  is,  an  arrangement  of  piezometric 
devices  might  be  able  to  provide  an  estimate  of  the  location  of  the  phreatic 
surface  over  a  period  of  time.  In  the  long  run,  nonintrusive  measurement 
procedures  that  can  detect  moisture  migration  would  be  highly  desirable. 

49.  An  initial  goal  of  physical  testing  for  degree  of  saturation 
should  be  the  validation  of  current  theories  and  models  for  moisture  migra¬ 
tion.  This  goal  can  probably  be  approached  through  small-scale  or  laboratory 
simulations  that  consider  both  seepage  and  capillary  action.*  Validation  of 
an  existing  or  new  model  for  moisture  migration  would  be  an  important  step  in 
improvement  of  service  life  in  concrete. 


*  or  diffusion;  see  e.g.  Bremer,  H.  W.  1965.  "Moisture  Migration  -  Concrete 
Slab -on -Ground  Construction,"  Journal  of  the  PCA  Research  and  Development 
Laboratories,  May,  pp  2-17,  also  PCA  Development  Dept.  Bulletin  D  89  and 
references  given  therein;  also  Powers,  T.  C.;  L.  E.  Copeland,  and 
H.  M.  Mann.  1959.  "Capillary  Continuity  or  Discontinuity  in  Cement 
Pastes,"  ibid  May  pp  38-48,  also  PCA  Research  Laboratory  Bulletin  110  and 
references  given  therein. 
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PART  VII:  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Summary 

50.  A  procedure  for  predicting  the  service  life  of  nonair-entrained 
concrete  subject  to  damage  due  to  freezing  and  thawing  has  been  developed 
which  addresses  both  the  known  and  uncertain  qualities  of  the  relevant  mate¬ 
rial  properties,  environmental  factors,  and  model  of  degradation  due  to  freez¬ 
ing  and  thawing  using  a  probabilistic  method.  Two  important  characteristics 
of  this  procedure  are:  (1)  it  rationally  addresses  the  uncertainties  inherent 
in  degradation  of  mass  concrete  due  to  freezing  and  thawing,  and  (2)  the  pro¬ 
cedure  is  mathematically  straightforward  for  implementation  by  CE  offices. 

51.  This  procedure  was  demonstrated  by  hindcast  application  to  two 
structural  features  at  CE  Civil  Works  structures  which  exhibited  an  apprecia¬ 
ble  degree  of  measurable  damage  due  to  freezing  and  thawing.  Required  data 
for  application  of  the  procedure,  e.g.  temperature  and  concrete  properties, 
were  available  for  these  features,  which  were  representative  of  typical  CE 
projects.  From  six  candidate  structures,  two  features  were  selected  as  case 
studies  for  application  of  the  procedure.  These  two  features  were  the  middle 
wall  and  land  wall  at  Dashields  Lock. 

52.  The  procedure  was  applied  using  the  available  data  for  each  case 
study.  One-  and  two-dimensional  thermal  analyses  were  used  to  determine  the 
thermal  response  of  each  wall  and  the  resulting  probabilities  of  critical 
temperature.  Simplified  seepage  analyses  provided  the  probabilities  of  crit¬ 
ical  saturation.  The  annual  probability  of  damage  and  the  predicted  service 
life  were  determined  from  the  joint  probabilities  of  critical  temperature  and 
critical  saturation  throughout  the  structure. 

53.  Sensitivity  of  concrete  service  life  prediction  to  the  uncertain¬ 
ties  of  input  factors  was  discussed.  Implications  of  physical  test  proce¬ 
dures,  and  discussion  of  future  needs  for  concrete  service  life  prediction 
were  examined. 
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Conclusions 


54.  Several  important  conclusions  can  be  reached  from  this  study. 
Current  procedures  for  thermal  modeling  and  analysis  appear  quite  adequate 
for  predicting  temperatures  in  a  concrete  structure.  Although  2-D  analyses 
are  better  for  determining  complex  thermal  response,  in  many  cases  a  series 
of  much  simpler  1-D  analyses  provide  a  very  good  estimation  of  temperatures. 
The  external  temperature  inputs  to  a  thermal  analysis,  i.e.  water  and  air 
temperatures,  were  shown  to  be  well  represented  by  sinusoidal  curves.  In 
this  respect,  the  internal  thermal  response  of  a  unit  volume  of  concrete 
likewise  follows  a  sinusoidal  variation.  This  observation  permits  a  simpli¬ 
fied  representation  of  the  inputs  and  results  of  the  thermal  response,  and 
also  simplifier  application  of  1-D  analyses.  An  examination  of  the  tempera¬ 
ture  distribution  for  any  point  in  the  structure  indicated  that  the  primary 
contributor  to  uncertainty  of  response  was  the  annual  variation  of  tempera¬ 
ture.  This  important  observation  permits  the  probabilistic  analysis  to  use 
the  results  of  only  a  single  thermal  analysis  at  the  mean  values  of  the  input 
factors . 

55.  The  general  understanding  and  analytical  models  for  predicting 
moisture  migration  and  degree  of  saturation  are  not  as  well  developed  as  for 
the  thermal  problem.  A  seepage  model  for  predicting  the  degree  of  saturation 
appears  to  provide  adequate  answers  for  the  prediction  of  service  life;  how¬ 
ever,  further  study  is  appropriate  to  substantiate  this  indication. 

56.  The  procedure  developed  and  presented  herein  provides  excellent 
agreement  with  observed  damage  due  to  freezing  and  thawing  at  the  two  sites 
studied.  The  general  trends  of  location  and  spatial  variation  of  damage  are 
very  similar  to  observations  and  measurements  at  the  two  sites.  More  encour¬ 
agingly,  the  actual  magnitudes  of  damage  predicted  by  the  procedure  compare 
favorably  with  the  previous  measurements.  This  result  provides  the  strongest 
indication  that  the  procedure  is  rational  and  would  enhance  the  ability  of  the 
CE  to  predict  service  life  at  its  many  other  concrete  structures. 

Recommendations 

57.  This  study  provides  a  strong  basis  for  prediction  of  concrete 
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service  life  at  typical  CE  structures,  but  would  be  strengthened  considerably 
by  investigating  several  areas  exposed  by  the  study.  Based  on  the  extremely 
favorable  results  of  the  two  case  studies  of  this  study,  the  procedure  should 
be  applied  to  several  other  structural  features  typical  to  CE  structures  for 
which  actual  damage  is  reasonably  well  known.  These  additional  applications 
of  the  procedure  would  likely  increase  the  confidence  in  the  procedure  and 
indicate  any  areas  that  require  refinement  before  broad  application  to  struc¬ 
tures  without  damage  measurements.  As  a  precursor  to  these  additional  appli¬ 
cations,  a  small  survey  of  structural  features  and  extent  of  damage  due  to 
freezing  and  thawing  would  indicate  the  most  appropriate  structures  to 
consider . 

58.  The  current  understanding  (or  lack  of  understanding)  of  moisture 
migration  in  concrete  structures  indicates  some  further  development  of  models 
for  moisture  migration  is  warranted.  This  development  would  include  more 
extensive  review  and  examination  of  the  current  theories  and  their  applica¬ 
bility  to  the  problem  at  hand. 

59.  As  a  natural  adjunct  to  the  further  refinement  of  the  procedure, 
the  procedures  for  physical  measurements  of  temperature  and  moisture  should 
be  extended  to  provide  direct  information  regarding  structure  state  and  vali¬ 
dation  for  the  analytical  models.  This  extension  is  relatively  straightfor¬ 
ward  for  thermal  measurements  but  will  likely  require  some  innovative 
improvements  for  measurement  of  moisture. 
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Table  1 

Constants  From  Sinusoidal  Fit  of  Temperature  Data 


Linear  Regression  of  Air  Temperature  Data  (Equation  4) : 
Air  temperature,  9a  -  A  +  B  sin(wt)  +  C  cos  (wt) 

-  A  +  D  sin  (wt  +  <f> ) 

A  -  50.275°  F 

B  -  -7.972°  F 

C  -  -21.187°  F 
D  -  22.637°  F 

Linear  Regression  of  Water  Temperature  Data  (Equation  5) : 
Water  temperat  re,  0w  -  E  +  F  sin(wt)  +  G  cos  (wt) 

-  E  +  H  sin  (wt)  +  <f>) 

E  -  56.538°  F 

F  -  -9.387°  F 

G  -  -20.018°  F 
H  -  22.110°  F 


Table  2 

Horizontal  Depth  of  Freeze -Thaw  Damage  in  Dashields  Lock  Walls 


Wall 

River 


Middle 


Typical 

Lock  Wall 

Lower  Gate 

Recess 

Side  of 
Wall 

Monolith 

No. 

Depth  of 
Damage 
ft 

Monolith 

No. 

Depth  of 
Damage 
ft 

Land 

L- 18 

0.0 

L-33 

1.5 

L- 19 

0.2 

L-  34 

1.2 

L-24 

0.0 

L-35 

1.3 

L-28 

0.0 

L-  30 

??? 

River 

M- 11 

0.0 

M-22 

0.6 

M-12 

0.1 

M-23 

2.8 

M-16 

0.1 

M-24 

2.9 

M-19 

0.8 

Land 

M-12 

??? 

M-20 

0.3 

M-16 

??? 

M-21 

0.1 

River 

R-13 

1.2 

R-24 

0.4 

R-20 

0.6 

R-  24 

1.0 

Land 
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Figure  1.  Determine  mechanism  of  damage  by  freezing  and  thawing 


a.  Probability  density  function. 


b.  Reliability  function. 


Figure  2.  Distribution  of  material  service  life 
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Figure  4.  Saturation  distribution 


Figure  5.  Dashields  Lock  and  Dam 
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Figure  9.  Temperature  boundary  conditions  for  middle  wall 
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Periodic  temperature  variation 
on  two  opposite  surfaces 


X(x')  = 


cosh  kx’  M+i ) 
cosh  kL*  (1+i) 


k2  =  w/2h2 


0)  =  27C/365 


x'  =  x  -  U2 


L*  =  L/2 


temperature  attenuation  =  |  X(x)  |  = 


V 


cosh  k  (2x-L)  ±  cos  k  (2x-L) 
cosh  kL  +  cos  kL 


p 

where:  h  =  thermal  diffusivity 

0e  =  external  temperature 
to  =  temperature  frequency 
L  =  wall  thickness 


Figure  12. 


One -dimensional  thermal  analysis  model  of  wall 


Temperature  attenuation,  X(x) 


distance  from  face,  ft. 

Figure  13.  Attenuation  of  temperature  with  distance 
and  thermal  diffusivity 


center 


temperatures  plotted  every  2  feet 
from  surface  to  center  of  wall 


bb  H 


50  H 


45  H 


30  -I 


25 

Jan 


_T_ 

Feb 


n - 1 - 1 - r 

Mar  Apr  May  Jun 


-i - 1 - r 

Jul  Aug  Sep  Oct  Nov  Dec 


month 


Figure  14.  Annual  temperature  cycles  at  various  depths 

(above  waterline) 
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Figure  17.  Comparison  of  temperatures  from  1-  and  2-D  analyses 

(2  feet  above  waterline) 
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Figure  18.  Comparison  of  temperatures  from  1-  and  2-D  analyses 

(2  feet  below  waterline) 
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Figure  19.  Comparison  of  probability  of  freezing 
for  1-  and  2-D  analyses 
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Figure  20.  Area  of  middle  wall  studied  in  more  detail 
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Figure  21.  Mean  temperatures  (°F)  in  middle  wall. 
See  Figure  20  for  orientation  of  this  section 


Figure  22.  Mean  temperature  (°F)  amplitudes  in  middle  wall. 
See  Figure  20  for  orientation  of  this  section 
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Figure  23.  Probability  of  freezing  in  middle  wall. 
See  Figure  20  for  orientation  of  this  section 


$  *  total  fluid  head  or  potential 
-  z  +  p/7 
where: 

p  -  fluid  pressure 
Y  -  fluid  unit  weight 
z  -  elevation  above  datum 


Figure  24.  Seepage  model  of  middle  wall 
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Figure  25.  Distribution  of  lock  pool  elevations 
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Figure  26.  Probability  of  critical  saturation  in  middle  wall. 
See  Figure  20  for  orientation  of  this  section 
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Figure  28.  Service  life  (yrs)  prediction  for  middle  wall  at  Dashields  Lock 
See  Figure  20  for  orientation  of  this  section 


Figure  29.  Section  through  land  wall  at  Dashields  Lock 
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Figure  30.  Temperature  boundary  conditions  for  land  wall 
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Figure  32.  Thermal  finite-element  model  of  land  wall 
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Figure  33.  Mean  temperatures  (°F)  in  land  wall. 
See  Figure  30  for  orientation  of  this  section 
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Figure  35.  Probability  of  freezing  in  land  wall. 
See  Figure  30  for  orientation  of  this  section 
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Figure  36.  Water  table  on  land  side  of  land  wall 


Figure  37.  Distribution  of  phreatic  surface  in  land  wall 
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Figure  38.  Probability  of  critical  saturation  in  land  wall. 
See  Figure  30  for  orientation  of  this  section 
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Lgure  39.  Annual  probability  of  damage  for  land  wall  at  Dashields  Lock. 
See  Figure  30  for  orientation  of  this  section 
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Figure  40.  Service  life  (yrs)  prediction  for  land  wall  at  Dashields  Lock 
See  Figure  30  for  orientation  of  this  section 


APPENDIX  A:  SINUSOIDAL  CURVE  FITTING 


For  cyclic  data  with  "n"  data  points  evenly  spaced  over  a  period  "T"  and  an 
equation  of  the  form: 

X  (t)  -  A  +  B  sincot  +  C  costot  ,  £0  -  2«/T 
the  coefficients  A,  B,  and  C  can  be  determined  by  minimizing 

Z  [  (Xj_  -  (A  +  B  sincot j_  +  C  coscot^)]2 
The  linear  regression  equations  may  be  written  as 


Zxi 

Zx^sincati 

ZXiCO3(0ti 

Zsinoat^  -  Zcosoat^  -  Zcostotisinoat^  -  0 

Zsin2tot^  -  /sin2(0tdt  -  T/2 

Zcos2<uti  “  Jcos2oatdt  -  T/2 

the  coefficients  are  determined  as 

A  -  (ZXj_)  /n  -  mean 

B  -  (2/T)Zxi3infflti 

C  -  (2/T)Zxico3(0ti 

For  the  alternate  equation  of  the  form: 

X(t)  -  A  +  D  sin  (cot  +  <|» 

then  D  -  (B2  +  C2)1/2 

<)>  -  tan_1(B/C) 


Zsinoati 

Zcoscot^ 


Zs  inoat 

Zsin2oati 

Zcosojt^sinoot^ 


Zcosoot^ 
Zcoscotj_sinooti 
Zcos2o)t  j_ 


Noting  the  following  simplifications  for  this  case : 


A3 


APPENDIX  B: 


SAMPLE  CALCULATIONS  FOR  TEMPERATURE  DISTRIBUTION 


D«^2w>x^rc4,  cut  d^{^\  <fj^  (o  ^0>£  -pt")  J^v  o**\ (a u^wt  ou_v  ”t"®wi  yj  eva-^t<-*^o  , 

SuLW'-jx-oi,  fe^ve^<ujfu^e  $_  =  $£  e£u*^  -  A  +  D  *°*  Ui^ 

/A  *  5o.2'75"  O^J  -  £.7^*  A_=  -p- 

5=  22.637°  <r0-  0.3IH-'  w=  a%4S 

h*°  /. 2o  f1/^  <>v=  <5-2-0 

De^exorvu^vwi  fx_o  Xe^^jiv'c-J'M^erf  (a)  *  =  0.  5  -^". 

&(*)~  A  +  X(x)  D  Sl^  cot 


h^l.Z  ^  k.  =  0.  0^6  4  -  Z./2737  kL'---2.o 3?74 

^1,0,  k=0.on  77  kl-  2.22659  kL'  — 2-/33?/ 

/-.  |,4  ,  * -0,07*?  I  kL  =  I,**/?!  /cL'~-J.fo34 

Xfr-U)  -  JT^l  . 

Y  (2.l2l37jr  c«  (Z.)ZniT)  1  3.7277 


X^M.o)  = 


e«4,  (-2.ns8l)  +c<ki  (-2.13381) 
c^Cv(2.22659)i-c^  (2.226S9) 


S  0,  9&  1? 


X(uhmV  j^t^zseet  . 

Y  C<h4v  (].??!?  I )  *C*-  (l.is/gf) 


B3 


S/JttPig  QL<u4un~toM 


Tou^Ur  S^rt^ 


9,(<)=  0„6>Aft)  =  50.  ns' 

46)  •  0o(M,6)*  X<«) -D  --  <3.962^6^2.637)  -  2/.7?  ' 

0;!6)  '  V*  0.765  2  =  0.570* 


4 0 f  ra)  -  X6)6*td)  *  0.762*4  (^s)-  ZZ.oli  ” 

4/ &-«•„)  -  x6)(b-<r0)  =  0-9bZ(,/,(U.ZZl)  =  2/.481  0 

4<T*v)  =  xi!=w)  5  =  .  ‘H.S*t>(n.bYT)*  2,.KH  • 
^ (p,!. ,  ' \ ,  x6‘- '•»)  0  ’  .Ksng(zz.(,n)  =  2/.707" 


■  -a^-y ,  /iiftVfttty' 


O.  0*176, 


Tev^p«^tU«  Dii+viluJb^v,  ; 

5  *  6  M*  50.  27S  “ 

o  ^ 

V’  <r02*  i  (V  *- 4J) 

=  .54  +  2  (.0916.  +  21-79  /  »  Z3?.oV  * 

<V  (5.43  " 


Table  B-l 

One  dimensional  temperature  distribution  above  waterline 


dist. 

annual  mean  Var(amplitude) 

mean  amplitude 

Variance 

total 

std.  dev. 

StA2 

divided 

X 

SoA2 

■5*SaA2 

•5*TaA2 

StA2 

21 

bv  total 

0.00 

0.59 

0.05 

256.22 

256.86 

16.03 

0.9975 

0.25 

0.59 

0.05 

246.52 

247.16 

15.72 

0.9974 

0.50 

0.59 

0.05 

237.44 

238.07 

15.43 

0.9974 

0.75 

0.59 

0.04 

228.95 

229.58 

15.15 

0.9973 

1.00 

0.59 

0.04 

221.01 

221.65 

14.89 

0.9971 

1.50 

0.59 

0.04 

206.73 

207.36 

14.40 

0.9970 

2.00 

0.59 

0.04 

194.38 

195.01 

13.96 

0.9968 

2.50 

0.59 

0.04 

183.77 

184.40 

13.58 

0.9966 

3.00 

0.59 

0.04 

174.73 

175.36 

13.24 

0.9964 

4.00 

0.59 

0.04 

160.70 

161.32 

12.70 

0.9962 

5.00 

0.59 

0.04 

151.07 

151.70 

12.32 

0.9958 

6.00 

0.59 

0.04 

144.80 

145.43 

12.06 

0.9957 

7.00 

0.59 

0.04 

141.00 

141.63 

11.90 

0.9956 

8.00 

0.59 

0.04 

138.92 

139.54 

11.81 

0.9956 

9.00 

0.59 

0.04 

137.93 

138.55 

11.77 

0.9955 

10.00 

0.59 

0.04 

137.56 

138.18 

11.76 

0.9955 

11.00 

0.59 

0.04 

137.47 

138.10 

11.75 

0.9954 

12.00 

0.59 

0.04 

137.47 

138.09 

11.75 

0.9955 

Table  B-2 

One  dimena ional  temperature  distribution  below  waterline 


dist. 

annual  mean  Var(amplitude) 

mean  amplitude 

Variance 

total 

std.  dev. 

StA2 

divided 

S 

SoA2 

.5*TaA2 

31£ 

£1 

by  total 

0.00 

3.87 

0.01 

244.42 

248.30 

15.76 

0.9844 

0.25 

3.87 

0.01 

235.17 

239.05 

15.46 

0.9838 

0.50 

3.87 

0.01 

226.51 

230.38 

15.18 

0.9832 

0.75 

3.87 

0.01 

218.40 

222.28 

14.91 

0.9825 

1.00 

3.87 

0.01 

210.84 

214.72 

14.65 

0.9819 

1.50 

3.87 

0.01 

197.21 

201.09 

14.18 

0.9807 

2.00 

3.87 

0.01 

185.43 

189.31 

13.76 

0.9795 

2.50 

3.87 

0.01 

175.31 

179.19 

13.39 

0.9783 

3.00 

3.87 

0.01 

166.68 

170.56 

13.06 

0.9773 

4.00 

3.87 

0.01 

153.30 

157.18 

12.54 

0.9753 

5.00 

3.87 

0.01 

144.11 

148.00 

12.17 

0.9737 

6.00 

3.87 

0.01 

138.14 

142.02 

11.92 

0.9727 

7.00 

3.87 

0.01 

134.51 

138.39 

11.76 

0.9720 

8.00 

3.87 

0.01 

132.52 

136.40 

11.63 

0.9716 

9.00 

3.87 

0.01 

131.57 

135.46 

11.64 

0.9713 

10.00 

3.87 

0.01 

131.22 

135.11 

11.62 

0.9712 

11.00 

3.87 

0.01 

131.14 

135.03 

11.62 

0.9712 

12.00 

3.87 

0.01 

131.14 

135.02 

11.62 

0.9713 
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